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Nonlocal effect is of fundamental importance to the un- 
derstanding of superconductivity, as it plays a key role 
in various physical phenomena observed in superconduc- 
tors. In general, the term "nonlocal" bears a meaning 
that the electromagnetic response of superconductors is 
not described by local equations due to the fact that 
the Cooper pairs carrying the supercurrent have a fi- 
nite distance of correlation (the "coherence length", ^o)- 
Meanwhile, the same term often refers to the situation 
that the physical quantities controlling superconductiv- 
ity depends on electronic momenta. In this short note, 
we discuss the nonlocal effect in the former sense, where 
it is pointed out that the effective magnetic penetration 
depth (A) in the multigapped superconductors may ex- 
hibit decrease with decreasing external field (i?o) in a 
manner A oc \/ Bq/ B*2 (with B*2 being the upper criti- 
cal field for the smallest energy gap) due to the nonlocal 
effect. 

In type I superconductors, it is well established that 
the nonlocal effect leads to the enhancement of the pen- 
etration depth, 

A-(A!Co)^/^ (1) 

where Al corresponds to the local London penetration 
depth. ^ The above estimation is obtained from the ele- 
mentary argument for a bulk superconductor, consider- 
ing a mean vector potential (A = XBq) over a surface 
layer of thickness A (which is zero elsewhere). Here, let 
us follow the similar line to clarify the situation in type 
II superconductors. We introduce a vector potential for 
an isolated flux line. 
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where $o is the flux quantum and r is the distance from 
the center of the flux. The above equation may be ap- 
proximated into a mean potential. 
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within a cylinder (0 r < A) and zero elsewhere (with 
S being the unit surface of flux line lattice). The average 
current density J around the flux can be estimated by 
using the nonlocal form of the supercurrent and A, 



where the factor A^/^q comes from the ratio of volume 
for the region of non-zero A (oc ttA^) to the effective 
integration volume (oc tt^q)- Combining eq. (4) with a 
Maxwell equation around the flux, we have 



Bo 1/ 4^ J A SBo 

-..|(curlB)| = — 

from which we obtain a solution 
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Using an expression for the upper critical fleld 
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eq. (6) is rewritten as 



A ~ At 
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Thus, when the nonlocal effect is significant (i.e., Al < 
^o), the penetration depth around flux line is renormal- 
ized by a factor \/Bo/ Bc2 ■ Interestingly, the effect dis- 
cussed here serves to reduce the penetration depth, which 
is in marked contrast to the case of eq. (1). 

It would be needless to mention that the situation con- 
sidered in the above estimation cannot be realized in or- 
dinary type II superconductors with a single gap, because 
the condition of Al < is never satisfied. However, in 
the case of multigapped supercondutors, it might happen 
that the coherence length corresponding to the smaller 
gap may be comparable to the penetration depth. To 
examine the practical situation, let us discuss an order 
parameter having two energy gaps, {i = a, tt) with 
Att < Act, that is actually realized in MgB2.^ The pres- 
ence of two gaps implies that there are two corresponding 
coherence lengths and associated upper critical fields 
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(i = CT,7r), 



(9) 
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with v-p being the Fermi velocity. Provided that Al < Ctt , 
the nonlocal effect leads to the renormalization of the 
penetration depth into an effective one. 
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A ~ At 
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where 6^ = 3^/3^2(7^)- It must be noted that super- 
conductivity is maintained by the larger gap (Aq-) for 
Bq < Bc2(a)i while A^r collapses and associated non- 



local effect disappears when Bq > B, 



c2(7r) 



Meanwhile, 



eq. (11) would not be valid for \ < ^a, where J is re- 
duced to zero. Thus, the nonlocal effect due to the smaller 

gap is predicted to occur only over the lower field range 

e^c2(.)/AE < Bo < i?c2(.) (or e/AE <K<1)- 

It is well established that muon spin rotation (^SR) 
measurement in type II superconductors provides a di- 
rect information on the spatial distribution of magnetic 
field 



B(r) = ^6(K)exp(-zKT), 



(12) 
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where K arc the vortex reciprocal lattice vectors. The 
Fourier component of the field profile, 6(K), is deter- 
mined by the penetration depth and coherence length; if 
we adopt the London model, we have 



6(K) = 



B^eM-UlK^) 



(13) 



1 + X^K"^ 

where (oc ^g) is the cutoff parameter to correct the 
nonlocal electrodynamics near the vortex cores. In a 
Gaussian approximation, the spin relaxation observed by 
/xSR is described by the Gaussian decay with a relaxation 
rate 

G{h) 
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(6 = Bo/B,2) (14) 



in which 7^ is the muon gyromagnetic ratio, and the 
factor G{h) cx (1 - 5)[1 + 3.9(1 - hf^/^ represents the 
reduction of mainly due to the stronger overlap of 
vortices (oc 1 — b) and additional decrease due to the 
contribution of vortex cores at higher fields.'^ The non- 
local effect for the two energy gaps may be incorporated 
by substituting Eq. (6) to the above equations so that 
the relaxation rate exhibits a field dependence 



\G(6)/A2, K 6x < Bc2(a)/i?c2W 



(15) 



where h = Bo/Bc2(a)- It implies that may exhibit 
strong deviation from that for the single gap below 
Bc2{tt) — *I'o/(27r^^) with a sharp increase with decreas- 
ing field (cr^ oc l/6,r). 

The predicted behavior of is surprisingly close to 
that found in the earlier reports on MgB2.''^"^ Prior to 
the establishment of double gap superconductivity, the 
strong enhancement of cr^ at lower fields was attributed 
to uncontrolled influence of flux pinning.^' ^ Although the 
flux pinning still remains as a possible cause of enhanced 
C7ju, we argue that the nonlocal effect discussed above 
may play a significant contribution to the observed field 
dependence of A. To visualize the situation, we repro- 
duce the previous data^ and a result of curve fitting by 
Eq. (15) with i?c2(7r) as a free parameter; the data below 
0.1 T were excluded from the fit as the field dependence 
of CTju was rciversccl probably due to the reduced flux den- 
sity near the lower critical field. The curve is in excellent 
agreement with data with B^2(Tr) — 0.53(1) T [^^ — 25(1) 
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Fig. 1. Muon spin relaxation rate in the superconducting state 
of MgB2 deduced by fits using a Gaussian damping (data re- 
produced from Ref.^). Sohd curve is a fit by the model taking 

account of nonlocal effect arising from double gap structure in 
the order parameter (see text for detail). 



nm] , reproducing the tendency of steep upturn with de- 
creasing field below the presumed upper critical field for 
the smaller gap. The field dependence of reported by 
other groups exhibit qualitatively similar trend with a 
sharp upturn of cr^ below ~0.5 T, and thereby it is likely 
that the relevant feature stems from an intrinsic property 
that might be ascribed to the nonlocal effect. 

The present discussion may cast some doubts to the 
earlier model^ for explaining the behavior of cr^ in MgB2 ; 
they assume that B{r) consists of a weighted sum of the 
London expression [Eq. (13)] with the cutoff replaced by 

and for the respective; components. Unfortunately, 
the London approximation would certainly fail for the 
smaller gap over the field range 6,r > 0.25.^ 

It must be noted that, while the agreement between 
data and the present model in Fig. 1 is encouraging, 
our estimation is no more than a very crude approxi- 
mation [e.g., the use of A in Eq.(4)] and more precise 
evaluation is clearly needed. Moreover, there are other 
sources of nonlocal effects that should be considered as 
well. For example, recent calculation based on the qua- 
siclassical Eilenberger theory indicates the presence of a 
nonlocal effect even for the single gap, where the effect 
may be observed as a gradual increase of A with increas- 
ing field. ^ A similar model for the multigap superconduc- 
tors is strongly awaited for the detailed understanding of 
the nonlocal effect. 
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